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Abstract 

Integrated engine-airframe optimal control technol- 
ogy ma y significantly improve aircraft performance. 
This technology requires a reliable and accurate pa- 
rameter estimator to predict unmeasured variables. To 
develop this technology base, NASA Dryden Flight 
Research Facility (Edwards, CA), McDonnell Aircraft 
Company (St. Louis, MO), and Pratt & Whitney (West 
Palm Beach, FL) have developed and flight-tested an 
adaptive performance seeking control system which op- 
timizes the quasi-steady-state performance of the F-15 
propulsion system. This paper presents flight and 
ground test evaluations of the propulsion system pa- 
rameter estimation process used by the performance 
seeking control system. The estimator consists of a 
compact propulsion system model and an extended 
Kalman filter. The extended Kalman filter estimates 
five engine component deviation parameters from mea- 
sured inputs. The compact model uses measurements 
and Kalman-filter estimates as inputs to predict un- 
measured propulsion parameters such as net propulsive 
force and fan stall margin. The ability to track trends 
and estimate absolute values of propulsion system pa- 
rameters was demonstrated. For example, thrust stand 
results show a good correlation, especially in trends, 
between the performance seeking control estimated and 
measured thrust. 

Nomenclature 

A, B, C, state variable model matrices 
A L, M 

AAHT high-pressure turbine area component 
deviation parameter, in 2 

* Aerospace Engineer. 

** Aerospace Engineer. Member AIAA. 

Copyright ©1992 by the American Institute of Aeronau- 
tics and Astronautics, Inc. No copyright is asserted in the 
United States under Title 17, U.S. Code. The U.S. Govern- 
ment has a royalty-free license to exercise all rights under 
the copyright claimed herein for Governmental purposes. All 
other rights are reserved by the copyright owner. 


AJ 

nozzle throat area, in 2 

BLD 

bleed air, Ib/sec 

CIVV 

compressor inlet variable guide vane 
angle, deg 

CPSM 

compact propulsion system model 

DEEC 

digital electronic engine control 

DEHPT 

high-pressure turbine component 
deviation parameter, percent 

DELPT 

low-pressure turbine component 
deviation parameter, percent 

DNOZ 

nozzle drag, lbf 

DRAM 

ram drag, lbf 

DWFNA 

fan airflow component deviation 
parameter, lb/sec 

DWHPC 

high-pressure compressor airflow compo- 
nent deviation parameter, lb/sec 

EPR 

engine pressure ratio, PT6/PT2 

F 

steady-state variable model sensitivity 
matrix 

FG 

gross thrust, lbf 

FN 

net thrust, lbf 

FNP 

net propulsive force, lbf 

FTIT 

fan turbine inlet temperature, °R 

h 

pressure altitude, ft 

HPX 

power extraction, hp 

K 

steady-state Kalman gain matrix 

KF 

Kalman filter 

M 

Mach number 

MIL 

military 

N 1 

fan rotor speed, rpm 

N1C2 

fan rotor speed, corrected to station 2, rpm 

N2 

compressor rotor speed, rpm 


OPR 

engine operating pressure ratio, PT4/PT2 

P amb 

ambient pressure, lb/in 2 

PB 

burner pressure, lb/in 2 

PLA 

power lever angle, deg 

PS 2 

static pressure at engine face, lb /in 2 

PSC 

performance seeking control 

PT 

total pressure, lb/in 2 (used with suffixes; 
list follows) 

RCVV 

rear compressor variable guide vane angle, 
deg 

RI 

Reynolds index, 8/6 124 

SMF 

fan stall margin, percent 

SMHC 

high-pressure compressor stall margin, 
percent 

SSVM 

steady-state variable model 

SVM 

state variable model 

TMT 

composite metal temperature, °R 

TSFC 

thrust specific fuel consumption, sec -1 

TT 

total temperature, °R (used with suffixes; 
fist follows) 

u 

vector of control variables in the state 
variable model 

WACC 

digital electronic engine control calculated 
airflow, lb /sec 

WCFAN 

fan airflow, lb /sec 

WCHPC 

high-pressure compressor airflow, lb/ sec 

WF 

gas generator fuel flow, lb/ hr 

X 

vector of state variables in the state 
variable model 

y 

vector of output variables in the 
steady-state variable model 

a 

angle of attack, deg 

0 

angle of sideslip, deg 

6 

reference pressure ratio, ambient pressure 
to sea-level pressure 

c 

vector of component deviation parameters 
in the compact propulsion system model 

8 

reference temperature ratio, ambient 
temperature to sea-level temperature 


Subscripts 

b predicted trim values 

c corrected 

m measured 


Superscripts 


T transpose 

estimated value of variable 
' augmented Kalman filter variables 

Prefix 


A perturbation 

Suffixes, PW1128 engine station numbers, 
Fig. 2 


2 

2.5 

3 

4 

4.5 
6 

7 


fan inlet 
compressor inlet 
compressor discharge 
high-pressure turbine inlet 
low-pressure turbine inlet 
afterburner discharge inlet 
nozzle throat discharge 


Introduction 

Digital engine control and optimal control algo- 
rithms enable significant performance improvements of 
the integrated aircraft-propulsion system. Developing 
and applying this technology will contribute to both 
commercial and military applications by maximizing 
thrust and fuel efficiency and extending engine life. 
Most benefits are directly attributable to advances in 
real-time parameter estimation techniques. Such tech- 
niques enable the control system to recover latent per- 
formance from the propulsion system which, until now, 
has been unrealized. 

To develop this optimal performance technology, 
NASA Dryden, McDonnell Aircraft Company, and 
Pratt & Whitney have developed and flight-tested an 
adaptive performance seeking control (PSC) system. 1 
The objective was to optimize the quasi-steady-state 
performance of the F-15 propulsion system. The 
adaptive features of the PSC system are provided by 
the propulsion system parameter estimation routine, 
which automatically adjusts the onboard model to 
more closely match the engine hardware. 

The PSC system was developed with the follow- 
ing optimization modes: minimum fuel at constant 
thrust, maximum thrust, and minimum fan turbine 
inlet temperature (FT IT) at constant thrust. Sub- 
sonic flight testing of the PSC algorithm was con- 
ducted at NASA Dryden covering all three modes at 
part- and military- (MIL) power conditions. Flight re- 
sults indicate that substantial benefits were obtained 
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from the PSC algorithm, up to 15-percent increases in 
thrust, up to 100 °R reductions in turbine tempera- 
ture, and between 1- and 2-percent savings in specific 
fuel consumption. 2 These results rely upon an accurate 
engine parameter estimation. 

A preliminary PSC evaluation of the propulsion sys- 
tem parameter estimation routine was conducted with 
limited flight data and determined to produce “reason- 
able estimates.” 3 In this earlier investigation only post- 
flight models were analyzed; the effects of measurement 
biases, flight condition, and engine degradation were 
not evaluated. In a recently completed subsonic flight 
program, test points were specifically designed to mea- 
sure the effectiveness of a propulsion system parameter 
estimator. This paper presents subsonic flight test and 
thrust stand evaluations of the parameter estimation 
process for nonafterburner power settings. Represen- 
tative data axe analyzed for the effects of engine degra- 
dation, flight condition, and measurement biases on the 
estimator. 

Aircraft Engine Description 

The PSC program has been implemented on the 
NASA F-15 research airplane (Fig. 1), which is a modi- 
fied high-performance aircraft capable of speeds greater 
than Mach 2. Two PW1128 afterburning turbofan en- 
gines power the F-15. The aircraft has been modified 
with a digital electronic flight control system. Addi- 
tional information on the F-15 can be found in Ref. 4. 

The PW1128 engine used in this study is a low- 
bypass ratio, twin spool, afterburning turbofan tech- 
nology demonstrator derived from the F100-PW-100 
engine. The engine uses a full-authority digital elec- 
tronic engine control system (DEEC) that is similar to 
the current production F100-PW220 engine controller. 
The DEEC provides closed-loop feedback control of 
corrected fan speed ( N1C2 ) and engine pressure ra- 
tio (EPR) through the fuel flow (WF) and the noz- 
zle area (A J) respectively. The compressor inlet vari- 
able guide vane (Cl W) and rear compressor variable 
vane (. RCVV ) positions are scheduled on rotor speeds 
through open-loop control. The DEEC software has 
been modified to accommodate PSC trim commands, 
however, the normal DEEC control loops (i.e., NIC2 
and EPR) have not been modified. A more detailed de- 
scription of the PW1128 engine can be found in Ref. 5. 

Two PW1128 engines were evaluated during the sub- 
sonic phase of the program. Initial PSC testing was 
conducted with a recently refurbished engine (located 
on the right side of the aircraft). All algorithm veri- 
fication and early testing were done with this engine. 
All PSC design work was performed using models of a 
nominal engine, which represented the refurbished en- 
gine. Near the end of the subsonic flight test program, 


a very degraded engine was placed on the left side of 
the aircraft. The main degradation was in the high- 
pressure rotor; both the compressor and turbine had 
significant degradation, estimated to be approximately 
2 percent in both areas. 

Instrumentation 

A diagram of the PW1128 engine is shown in 
Fig. 2. The locations of the DEEC instrumentation, 
DEEC-calculated parameters, and the parameters e&- 
timated by the PSC system are indicated. Fan airflow 
(WCFAN) and engine face total pressure (PT2) were 
independently modeled by both the DEEC and PSC 
control laws. The PSC algorithm requires only con- 
ventional DEEC-instrumented parameters as inputs, 
and estimates other necessary parameters within the 
algorithm. The engine instrumentation and PSC pa- 
rameters were sampled at 20 samples/sec. 

In addition to the basic engine parameters detailed 
earlier, the challenging nature of the technology be- 
ing demonstrated required the recording of many in- 
ternal algorithm variables. These additional variables 
provided for real-time and postflight analysis or de- 
bugging of the algorithm. More than 200 internal pa- 
rameters were recorded at a rate of 100 samples/sec. 
In the Kalman-filter (KF) estimator, the inputs, out- 
puts, and residuals were recorded. At the compact 
propulsion system modeling stage, all estimated and 
measured inlet and engine parameters were recorded 
including instrumented temperatures, pressures, and 
control positions, and estimated stall margins, thrust, 
and drag components. In the optimization phase, 
the constraints, optimal solution, and optimizer health 
condition codes were recorded. Finally, the actual com- 
mands that were sent to the engine through the DEEC 
were recorded. 

The airdata used in the PSC were obtained from the 
F-15 production side probes. The algorithm corrects 
the data for position error and location effects. All data 
were recorded on a pulse-code-modulation system. 

Performance Seeking Control System 

The PSC algorithm is a control law that attempts to 
optimize propulsion system steady-state performance. 
The algorithm in Fig. 3 includes parameter estima- 
tion, modeling, and optimization. Flight measure- 
ments are used as input to the propulsion system pa- 
rameter identification process and then to the compact 
propulsion system model (CPSM), as described later. 
Once the updated engine model is obtained, a linear- 
programming optimization routine is used to generate 
trims to control parameters. An essential assumption 
of the estimation and modeling processes is the steady- 
state nature of the problem. According to program 
objectives, only existing engine control measurements 
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were used, even though performance could have been 
improved with additional measurements. 

An extended KF is used to identify the differences 
between the nominal model and the actual flight article 
by producing five engine efficiency parameters which 
represent these differences. The estimated parame- 
ters account for low- and high-pressure core efficiencies 
(DELPT and DEHPT , respectively), differences in 
fan and high-pressure compressor airflow (DWFNA 
and DWHPC respectively), and a change in high- 
pressure turbine area {AAHT). 

In the second stage of the estimation, flight mea- 
surements are augmented with the KF efficiencies to 
develop the CPSM. The CPSM uses the KF estimates 
to adjust the steady-state model to accurately repre- 
sent the flight article. 

Finally, the CPSM provides estimates of unmea- 
sured propulsion system outputs, such as net propul- 
sive force (FNP) and fan stall margin ( SMF ), to the 
follow-on optimization routine. The optimization rou- 
tine is based on the linear-programming simplex algo- 
rithm. For this paper, there were three optimization 
modes: maximum thrust, minimum FT IT at constant 
thrust, and minimum fuel consumption at constant 
thrust. The optimization sends the optimal trims to 
the DEEC, which applies trims to the engine. 

Kalman-Filter Estimator 

The adaptive feature of the PSC algorithm is pri- 
marily provided by a KF, which estimates engine effi- 
ciency. The KF estimates five engine efficiency parame- 
ters, which account for engine- to-engine variations and 
engine deterioration. Engine operating efficiency has 
been modeled in a state-space perturbation formula- 
tion as follows: 

Af = [A] Ax A [B]Au + [L] C 

Ay = [C] Ax + [D]Au + [M] C (1) 

where 

x = [N1 N2 TMT] t 

u = [WF AJ CIVV RCVV HPX BLD} T 

y = [PT6 PT4 FTIT N 1 N2] T 

C = [DELPT DEHPT AAHT DWHPC DWFNA] T 

The state, control, measurement, and efficiency vec- 
tors are x, u, y, and £ respectively. Horsepower ex- 
traction ( HP A) and bleed airflow i^BLD) are modeled 
as part of the control vector, and aggregate turbine 
metal temperature ( TMT ) is modeled as part of the 
state vector. The state, control, and measurement per- 
turbations are calculated from the engine data (x, u, 
y) and the predicted state, control, and measurement 
steady-state trim values (x*,, u b , y b ) where 


Ax = x — Xb 
Au = u — Ub 
Ay = y~y b 

In the PSC problem, the values for £ are unknown and 
require estimation. 

Engine deterioration occurs slowly relative to the dy- 
namics of the engine state vector and is assumed to be 
steady state. Hence, ( can be approximated to be 0. 
Reformulating the problem by augmenting the state 
vector with ( yields: 


Ax 


A L 1 

Ax 


" B ‘ 

. . 


.0 oj 

. A c . 

+ 

0 


[Au] 


(2) 


and 


[Ay] = [C M] 


Ax 

AC 


+ [D] [Au] 


After combining the Ax and ( vectors to form an aug- 
mented state vector, the above equations are written 
in a compact KF format yielding 


Ax' = [A'] Ax' + [B']Au + [K](Ay - Ay) 
Ay = [C’} Ax' + [D]Au 


where the KF gain matrix, K, is a function of the solu- 
tion to the steady-state matrix Ricatti equation. More 
details concerning the KF design may be found in Lup- 
pold et al. 6 

The block diagram in Fig. 4 shows how the KF 
estimator employs the dynamic state variable model 
(SVM). The SVM is a piecewise linear model en- 
compassing the entire range of engine operation at 
Mach 0.9 at an altitude of 30,000 ft at standard day 
conditions. It consists of a state-space perturbation 
model and an associated table of steady-state trim val- 
ues for all engine variables in the model. The pertur- 
bations of Ax, Ay, and Au represent the differences 
between flight-measured and model steady-state values 
(also referred to as basepoints, hence the subscript 6). 
The predicted steady-state trim values are stored as a 
bivariate function of PT4 and PT6 and are for a nomi- 
nal, undeteriorated engine. There are 49 sets of A' , B', 
C ', D , and K matrices corresponding to values of PT4 
ranging from 23 to 260 lb/in 2 , which accommodate the 
engine operating range for the flight envelope corrected 
to the 0.9-Mach, 30,000-ft altitude reference condition. 

Figure 2 shows the locations of the engine parame- 
ters. Values for the following measurements and con- 
trol variables are taken directly from flight data: Nl, 
N 2, burner pressure {PB), FTIT , PT6, WF, AJ 
CIVV, and RCVV. The PT4 is modeled as a func- 
tion of PB, HPX is modeled as a function of N2 , and 
BLD is modeled as a function of Mach and altitude. 
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Additional flight measurements are used by the KF for 
calculating other engine variables and to transform the 
engine data to the SVM design condition of Mach 0.9 
at an altitude of 30,000 ft. Standard correction factors 
are used to transform the measured and calculated en- 
gine variables to different flight conditions. 

The engine efficiency parameters convey operational 
differences between an engine model at 0.9 Mach and 
30,000 ft and the measured in-flight engine. The esti- 
mates were designed to represent the most significant 
influence on engine efficiency -namely, engine degrada- 
tion. Included in engine degradation are component 
deterioration and engine-to-engine variations. How- 
ever, from the mathematical formulation, the five esti- 
mates must also contain all other unknown effects, such 
as measurement bias and Reynolds effects. The source 
of measurement bias is sensor error. Reynolds ef- 
fects arise from differences between the design point at 
0.9 Mach and 30,000 ft and the flight-tested Mach and 
altitude. Effectively, the five estimates can be thought 
of as model-matching parameters or component devia- 
tion parameters. As such, they will be termed compo- 
nent deviation parameters throughout the remainder of 
this paper because they are more than simple operating 
efficiencies. 

Compact Propulsion System Model 

The second phase of the PSC parameter estimation 
algorithm formulates the CPSM. The CPSM is a sim- 
plified steady-state model of the engine and inlet. Es- 
timates from the KF are included in the control input 
vector to the CPSM to adjust the model to more accu- 
rately reflect engine operation. The CPSM consists of 
a linear steady-state perturbation model of the engine, 
steady-state trim tables, and follow-on nonlinear cal- 
culations. Figure 5 shows the structure of the CPSM. 

The linear portion of the CPSM is the steady-state 
variable model (SSVM). The SSVM has the form 

Ay = [P]Au (3) 

where 

Ay = 2/ - 2/6 

A u = u — Ub 

The symbols u and y represent the control input and 
measurement vectors, respectively. They are defined 
to be 

u = [WF PT6 CIVV RCVV HPX BLD DEHPT 
DELPT DWFNA DWHPC AAHT] t 
y = [N1C2 N2 AJ PT2.5 PT4 TT2.5 TT3 TT4 FTIT 
TT6 WCFAN WCHPC\ T 

Although this model was developed for 0.9 Mach at an 
altitude of 30,000 ft, it was translated to the sea-level 


static standard day reference condition using standard 
correction factors. The steady-state trim values are 
analogous to those used with the SVM in the KF, pro- 
viding the y b and u b vectors in equation (3) . The SSVM 
trim and the matrix models ( y b , u b . and F) are sched- 
uled as a bivariate function of PT4 and PT6 using 
linear interpolation between model points. 

The SSVM uses engine measurements for the follow- 
ing control inputs: WF, PT6, CIVV , and RCVV. 
The HPX and BLD are modeled as in the KF. The 
KF estimates are input to the SSVM calculation as 
part of the control vector. The SSVM provides esti- 
mates for the measurement vector y at the sea-level 
static condition and then transforms those estimates 
back to the actual flight condition with the standard 
correction factors. 

Following completion of the linear SSVM calculation, 
the nonlinear CPSM estimates are calculated. Refer- 
ring again to Fig. 5, these variables include PT7, TT7, 
gross thrust ( FG ), FNP, ram drag (DRAM), nozzle 
drag ( DNOZ ), AJNL, SMF , and high-pressure com- 
pressor stall margin ( SMHC ). The nonlinear calcu- 
lations use a combination of analytical equations and 
empirically derived data tables. They are based on 
measured engine variables and SSVM estimates. If 
a variable is both measured and estimated, the flight 
measurement is used in the nonlinear calculations. The 
nonlinear outputs from the CPSM are used with the 
SSVM outputs by the follow-on optimization routine. 

Performance Seeking Controls Optimization 

The PSC optimization routine is based on linear- 
programming techniques. As input the optimization 
uses a propulsion system matrix derived from the es- 
timated engine variables by the CPSM. The solution 
to the linear-programming problem represents the local 
minimum within the defined constraints. By definition, 
the optimal solution will always be on the boundary of 
at least two constraints. Depending on the optimiza- 
tion mode, either equality or inequality constraints are 
placed on parameters such as FTIT, FNP, and phys- 
ical boundaries. 

Of the possible 18 optimization constraints, the 
SMF estimate is often one of the two active con- 
straints. The optimization minimum SMF limit is 
determined by the accuracy of its calculation. The 
inputs to the SMF calculation are EPR and fan air- 
flow (WACC). Since EPR is considered very accu- 
rate, the WACC estimate drives the hmit for mini- 
mum SMF. Without the adaptive features of the PSC 
system, the error in the WACC estimate may be as 
high as ±4 percent leading to a required SMF hmit 
of approximately 10 percent. With the adaptive fea- 
tures of PSC, the increased accuracy in the WACC 
estimate leads to a lower required minimum SMF of 
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4 percent. With, the reduced constraint, the engine can 
be uptrimmed further. Therefore, the increased accu- 
racy of the adaptive parameter estimation process leads 
directly to improved engine and aircraft performance 
benefits. 

The PSC algorithm is closed-loop in a very general 
sense, since flight measurements are fed into the es- 
timation and modeling. However, the optimal solu- 
tion is based on models rather than flight data; as 
such, the algorithm may be described as a model-based 
command system. 

Evaluation Procedure 

The PSC subsonic flight testing and thrust stand 
testing provided an abundance of data with which to 
evaluate the propulsion system parameter estimation 
process. Maneuvers flown include aircraft accelerations 
and cruise flight at various altitudes and power settings 
on both the refurbished and degraded engines. Hence, 
a spectrum of flight conditions, engine operating condi- 
tions, and engine degradation levels was available from 
which to select data accurately representing the capa- 
bility of the PSC estimation process. 

The evaluation of the parameter estimation process 
involved analyzing the KF engine component deviation 
parameters first and then the CPSM estimated engine 
variables. The KF estimates have no associated mea- 
surement with which to compare. The interpretation of 
the absolute values of the estimates was somewhat lim- 
ited because the KF model is not a fully accurate rep- 
resentation of the actual propulsion system. The over- 
riding factors determining the value of the estimates 
and particularly their trends are power setting, engine 
degradation, measurement bias, and Reynolds effects. 
Of these known influences on the estimates, only power 
setting and engine degradation were modeled. Thus, 
some insight may be gained from observing trends in 
the estimates in comparison with engine degradation 
and power setting, or equivalently, operating pressure 
ratio (OPR, defined as the ratio of total burner pres- 
sure, PT4, to inlet total pressure, PT2, to normalize 
for power setting). However, since unmodeled sources 
of estimation error are not separable, it is not possi- 
ble to measure the sensitivity of each estimate to each 
known source of error. 

Evaluation of the CPSM parameter estimation is 
most readily accomplished by direct comparison with 
engine measurements. Several of the SSVM- estimated 
parameters are also part of standard FI 00 instrumen- 
tation: JV1C2, iV2, AJ , PT4, and FTIT. The param- 
eter jVlC2 is the measured fan rotor speed corrected 
to station 2 with TT2. In addition, the degraded test 
engine had special instrumentation for PT2.5, TT 2.5, 


and TT3. Comparisons were made between the DEEC 
airflow calculation and the CPSM-estimated airflow. 

Engine Testing 

Each of the three PSC optimization modes was 
tested on the refurbished and degraded engines. The 
minimum fuel mode was tested in steady- state cruise 
maneuvers, the minimum FTIT in steady-state cruise 
and acceleration maneuvers, and the maximum thrust 
mode in steady-state cruise and acceleration maneu- 
vers. Steady-state cruise test points usually consisted 
of 2 min with the PSC disengaged (PSC trims calcu- 
lated but not sent to the engine) followed immediately 
by 2 min with the PSC engaged (PSC trims calculated 
and sent to the engine). The 2-min time interval was 
determined to be an adequate period for the PSC al- 
gorithm to reach steady state. Typically, the KF es- 
timates required the most time to reach steady state, 
indicating a relatively large time constant in the filter. 
Before commencing the test, the power setting was ad- 
justed to the required position and was left there for the 
duration of the test maneuver. Acceleration test points 
required sustained turns for engine and algorithm sta- 
bilization after moving the throttle to the required test 
position. Once the engine-algorithm (i.e., the FTIT 
and PSC parameter estimates) stabilized, the pilot lev- 
eled the wings and began the acceleration test, d he 
baseline configuration tests were conducted at three 
altitude and Mach number conditions: 15,000 ft and 
0.9 Mach, 30,000 ft and 0.9 Mach, and 45,000 ft and 
0.88 Mach. 

Besides the testing completed in the nominal PSC 
configuration, some parametric testing was also per- 
formed to determine overall algorithm robustness. 
During the course of the subsonic flight testing, biases 
on engine measurements were adjusted several times to 
more accurately reflect the engine hardware and instru- 
mentation sensors. The effect of measurement biases 
on the estimation is of interest since it will demon- 
strate some of the sensitivity of the PSC algorithm to 
measurement errors. 

Degraded Engine Testing 

Tests were performed on the degraded engine to 
demonstrate the ability of the PSC system to adapt 
to various levels of engine deterioration. Specifically, 
the engine deterioration refers primarily to the high- 
pressure turbine. These blades and vanes were partially 
eroded from extended use (if the blades had been in a 
production engine they would have been replaced). To 
a lesser extent, the high-pressure compressor was also 
deteriorated, mainly because of worn tip seals. The 
poor condition of both components adversely affected 
engine performance. For example, the deterioration 
was manifested in FTIT ; the degraded engine operated 
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at hotter temperatures than the refurbished engine at 
the same flight condition and power setting. 

The state of the degraded engine is reflected in the 
KF engine component deviation parameters, CPSM en- 
gine estimates, and the optimal engine operating point. 
Assuming identical bias configuration and flight condi- 
tions, a comparison of the KF estimates between the 
refurbished and degraded engines will reflect relative 
differences in engine operating efficiency. Because of 
the deteriorated components, the state of the degraded 
engine differed from the nominal engine and resulted in 
different estimates and optimization results. Maneu- 
vers were conducted in the same manner as described 
for the refurbished engine testing. 

Thrust Stand Testing 

A thrust stand test was performed during the course 
of the subsonic flight testing to demonstrate and cali- 
brate real-time thrust calculations. 7 Degraded and re- 
furbished engines were tested on the aircraft with the 
PSC system engaged and disengaged. Effectively, the 
test setup consisted of the aircraft being tied down with 
a cable which was attached to strain gages to mea- 
sure force supplied by the engines. The strain gages 
measured the gross thrust as well as any wind-induced 
forces, or equivalently, net thrust. Although some 
thrust stand procedural or hardware error may have 
contributed to inaccuracy in the absolute thrust read- 
ings, the trends in thrust should be accurate. Thus, the 
thrust stand provides a unique opportunity to evaluate 
the onboard thrust estimate. 

Results and Discussion 

Results are presented to illustrate the ability of the 
parameter estimation process to provide acceptable es- 
timates of propulsion system outputs. Algorithm eval- 
uation of the overall estimation is approached at three 
levels: First, KF flight data will be qualitatively an- 
alyzed for trends and comparisons made with pre- 
dictions; second, the CPSM estimates will be com- 
pared with available measurements; third, the effect 
of parametric variations on the estimation process will 
be shown. 

Kalman Filter Component Deviation 
Parameters 

The trends and levels of the five KF estimates Eire 
evaluated at different flight conditions, power settings, 
and engine health. The values of DELPT , DEHPT, 
AAHT , DWHPC , and DWFNA are each plotted as 
a function of OPR. In general, the OPR points reflect 
tests done at nominal power level angle (PLA) values 
of 40°, 50°, 60°, 70°, and MIL. Engine component 
deviation parameters from refurbished engine testing 
at the KF design point of 0.9 Mach and 30,000 ft are 


presented to establish baseline magnitudes for these es- 
timates. Then, the variation of the engine component 
deviation parameters with respect to flight condition is 
evaluated. Next, the sensitivity of the KF estimation 
to measurement biases is examined. Finally, the effect 
of engine degradation on the KF estimates is presented. 

Besides the modeled engine degradation effects, 
the KF engine component deviation parameters in- 
clude many effects which are not explicitly modeled. 
Among the more influential factors are flight condition 
(Reynolds effects) and measurement bias. Boundary- 
layer and flow separation variations with flight condi- 
tion are examples of Reynolds effects. Measurement 
biases are known to exist for the engine instrumenta- 
tion of AJ, FT IT, and WF. Measurement bias esti- 
mates were obtained by the engine manufacturer based 
on detailed analysis of ground runs and flight results. 
Accordingly, the data presented in this paper were ob- 
tained with the resulting set of biases: +70° for FT IT, 
4-2.0 in 2 for AJ, and 4-180 lb/hr for WF for the re- 
furbished engine, and +70° for FT IT for the degraded 
engine. Thus, any remaining bias is probably small. 
Later, results are presented with a preliminary bias 
configuration (+40° for FT IT, -2.9 in 2 for AJ) used 
with the refurbished engine. 

As stated earlier, the PSC command system is model 
based, and as such, will produce similar results whether 
engaged or disengaged; differences are caused by rel- 
atively small changes in the engine operating point. 
Since the objective of this paper is to evaluate the es- 
timator, this evaluation will be comparable whether 
or not PSC is commanding optimal trims to the en- 
gine. The results obtained with the PSC system disen- 
gaged were more complete and of better quality; there- 
fore, these results are used for most results presented 
throughout the paper. 

Baseline Results 

Steady-state testing with the refurbished engine at 
0.9 Mach and 30,000 ft minimizes the influences of 
flight condition and engine degradation on the KF 
estimates. Still, there are modeling errors. The 
refurbished engine is not exactly the same as the en- 
gine model, and there are sensor biases. If the model 
contained no error, the KF estimates would be zero for 
the nominal engine at the design flight condition. 

Figure 6 presents the KF engine component devia- 
tion parameters for refurbished engine testing done at 
0.9 Mach and 30,000 ft. The OPR ranges from 12 at 
40° PLA to 29 at MIL PLA. Overall, the estimates 
are well behaved and follow smooth trends indicating 
reasonable KF estimates. However, not surprisingly, 
the onboard-determined engine component deviation 
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parameters are not zero. The following paragraphs dis- 
cuss the component deviation parameters of Fig. 6: 

Low-Pressure Turbine Component Deviation 
Parameter. The estimated DELPT (combination of 
fan and low-pressure turbine) is low relative to the en- 
gine model, ranging from —15 percent at OPR = 12 to 
— 1 percent at OPR — 29. The —15 percent at OPR = 
12, where PLA is 40°, may be too low since at other 
OPRs, the DELPT lies within ±5 percent. 

High-Pressure Turbine Component Deviation 
Parameter. The estimated DEHPT (combination 
of compressor and high-pressure turbine) of the refur- 
bished engine decreases very little relative to the en- 
gine model with increasing engine operating ratio, from 
0.2 percent at OPR — 12 to -3 percent at OPR = 29. 

High-Pressure Turbine Area Component De- 
viation Parameter. The estimated AAHT (or gas 
flowpath area) of the refurbished engine decreases rel- 
ative to the engine model with increasing engine oper- 
ating ratio, from 1.3 in 2 at OPR = 12 to —1.3 in 2 at 
OPR = 29. 

High-Pressure Compressor Airflow Compo- 
nent Deviation Parameter. The estimated 
DW 11 PC displays trends similar to those observed 
with the high-pressure spool component deviation 
DEHPT. The DW hi PC of the refurbished engine 
decreases very little relative to the engine model with 
increasing engine operating ratio, from 1 lb/sec at 
OPR = 12 to -1 Ib/scc at OPR = 29. 

Fan Airflow Component Deviation Parame- 
ter. The estimated DWFNA exhibits a minimum of 
-4 lb/sec at OPR = 18. The DWFNA of the refur- 
bished engine increases overall with respect to the en- 
gine model from —3 lb/sec at 40° PLA to —0.7 Ib/sec 
at MIL PLA. The 70° PLA point where OPR = 26 
appears to be another irregularity that is inconsistent 
with the general trend. 

An indication of the significance of the KF estimates 
is obtained by comparing the flight test and theoret- 
ically predicted values. From the data, the estimates 
were seen to have nonzero values indicating that the 
flight article is not like the model. Whether these dif- 
ferences are important w’ill be determined by examin- 
ing the effects of known biases and engine degradation 
at 0.9 Mach and 30,000 ft. To evaluate the impact 
of the magnitude of the KF engine component devia- 
tion parameters on the total PSC system, a compari- 
son is presented later for the CPSM estimation process 
with KF estimates as input and with zeroes as input 
at 0.9 Mach and 30,000 ft. 


Flight Condition Effects 

The KF estimation was evaluated with steady- 
state flight test data at three different altitude 
and Mach-number combinations for the degraded en- 
gine: 0.9 Mach, 15,000 ft; 0.9 Mach, 30,000 ft; and 
0.88 Mach, 45,000 ft. Corresponding to each flight con- 
dition is a Reynolds index (RI). The RI is calculated as 
a function of inlet stagnation temperature and pressure 
(0 = 7T2/518.9 and S = FT2/14.696; RI = <5/0 124 ), 
and is an indication of the thermodynamic properties 
of the free-stream air. The variation of each engine 
component deviation parameter with flight condition 
will indicate the overall effect of RI. Engine degra- 
dation and measurement bias are assumed not to vary 
with flight condition, but Reynolds effects certainly do. 
By examining each estimate at the three tested condi- 
tions, the Reynolds effects will be indicated. Each of 
the engine component deviation parameters is plotted 
as a function of OPR at the three flight conditions in 
Fig. 7. 

The parameter DELPT of the degraded engine 
shows a strong dependency on flight condition. For 
power settings ranging between an OPR of 10 and 30, 
the deviation is highest at 0.9 Mach and 15,000 ft and 
lowest at 0.88 Mach and 45,000 ft relative to the nom- 
inal engine. Differences from the 0.9-Mach, 30,000- 
ft design point are clearly seen in Fig. 7; at 15,000 
ft, DELPT ranges up to 7-percent higher, and at 
45,000 ft it varies as much as 7-percent lower. 

The estimated DWFNA and to a lesser extent 
DEHPT , like DELPT , also display tendencies with 
flight condition. Such tendencies are most noticeable 
at 15,000 ft and least noticeable at 45,000 ft. The esti- 
mated AAHT and DW 11 PC exhibit less pronounced 
trends with flight condition than the other three KF 
estimates. In general, AAHT and DW 11 PC show 
small increases with increasing altitude. Other than 
one test point at OPR = 24 at 45,000 ft, the estimated 
DW H PC is negative at all test conditions, indicating 
less high-pressure compressor airflow than is modeled 
by the nominal engine. The probable cause of the sub- 
nominal airflow 7 is engine degradation relative to the 
model. 

To determine the individual contributions of 
Reynolds effects to the estimation would require know r - 
ing the biases. The fact that these results were 
obtained with a degraded engine should not affect con- 
clusions relative to Reynolds effects. Of the three flight 
conditions tested, all five estimates displayed the least 
variation with OPR at the 0.9-Mach, 30,000-ft design 
condition. The DELPT , DEHPT , and DWFNA 
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were seen to decrease, and the A A FIT and DW II PC 
increase with increasing altitude or decreasing RI. 

Measurement Bias Results 

The refurbished engine was tested with two sets of 
measurement biases for AJ , WF, and FT IT . To cor- 
rect for what was believed to be measurement bias, 
corrections were added to the appropriate measure- 
ments before being used in the PSC algorithm. What 
will be referred to as the final bias set 1 consisted of 
+70° FT IT, +2.0 in 2 AJ , and +180 lb/hr WF mea- 
surement biases, and the preliminary bias set 2 con- 
sisted of +40° FT IT and -2.9 in 2 AJ. Data were 
collected at 0.9 Mach and 30,000 ft for cruise points at 
four partial- and MIRpower settings for the final bias 
set 1 and two partial- and MIL-power settings for the 
preliminary bias set 2. The resulting KF estimates are 
shown in Fig. 8 as a function of OPR. The data for bias 
set 1 represent the biases believed to be most accurate. 
Bias set 2 are values used early in the program. 

The estimated DELPT of the refurbished engine 
tends toward zero with increasing OPR for both bias 
sets. The two bias sets result in DELPT differing by 
up to 23 percent at the lowest OPR (~ 40° PLA) but 
by only 3 percent at the highest OPR (~ MIL PLA). 
For all tested OPRs ( 1 2~30) , the final bias set 1 re- 
sulted in lower values for DELPT. Similarly, the dif- 
ferences in the estimated DWFNA between the bias 
sets are greatest at the lowest OPRs (12 lb/sec) and 
least at the highest OPRs (1 lb/sec). 

The estimated DEH FT for the final bias set 1 re- 
sulted in higher values for DEHPT , over the range 
of OPRs. The differences in the estimated AAIIT 
between the bias sets decrease with increasing OPR. 
Overall, the final bias set 1 resulted in greater values 
for A A! FT . For all tested OPRs, bias set 1 resulted in 
more positive values for DW II PC. 

In general, it appears that the effect of biases is more 
exaggerated at lower OPRs or power settings and that 
the biases have less effect at the higher power settings. 
At MIL PLA, some engine parameters are running on 
operating limits and physical boundaries, and as such 
there is little or no effect of measurement bias. How- 
ever, as PLA is reduced, fewer engine parameters are 
operating on limits, and the effects of measurement 
bias will be more pronounced. KF estimation sensi- 
tivity to individual measurement bias on FT IT, AJ , 
or W F (as opposed to a combined bias set) is unclear 
since no parametric bias data are currently available. 
The results presented above represent the cumulative 
effects of a difference in biases of 30° FT IT, 4.9 in 2 AJ , 
and 180 lb/hr WF. Additional bias evaluations are 
made later for the CPSM estimation. 


Engine Degradation Results 

Steady-state data gathered on the degraded and re- 
furbished engines (with bias set 1) at 0.9 Mach and 
30,000 ft are compared to gauge the overall effect of 
engine deterioration on the KF estimates. The data 
are shown in Fig. 9. 

The estimated DELPT displays very little varia- 
tion between the refurbished and degraded engines 
except at lower OPRs. Estimates of DELPT for 
both engines increase with increasing OPR. The es- 
timated DEIIPT also shows small variations with en- 
gine degradation. Estimates of DEHPT for both en- 
gines decline slightly with increasing OPR. 

The estimated AAIIT decreases with increasing 
OPR for both engines, although there is some indi- 
cation of differences with engine degradation. The 
AAIIT estimate for the refurbished engine decreases 
by 3 in 2 over the OPR range, whereas the AAIIT esti- 
mate for the degraded engine decreases by only 0.5 in 2 . 

The estimated DW II PC decreases more for the re- 
furbished engine than for the degraded engine. For all 
tested OPRs, the degraded engine results show lower 
DW H PC estimates than those for the refurbished en- 
gine, as much as 3 lb/sec lower. The DWFNA vari- 
ation with OPR for the refurbished engine is much 
greater than that observed with the degraded engine. 

For the data in Fig. 9, the condition of the degraded 
engine is not obvious based on the engine component 
deviation estimates. This may indicate insufficient fi- 
delity in the KF model and suggests that the engine 
component deviation estimates are more accurately re- 
ferred to as model-matching parameters. Data clearly 
do not correlate with known degradation of the high- 
pressure spool. Comparisons of absolute magnitudes of 
any one KF estimate are inconclusive as to the effects 
of engine degradation, even at the design point of 0.9 
Mach and 30,000 ft. Overall, the KF estimation for the 
refurbished engine shows more sensitivity with respect 
to power setting than the degraded engine. Based upon 
the KF model, it would be expected that the trends in 
the engine component deviation parameters with OPR 
would be relatively constant or at least small. These 
results may indicate better overall modeling for the de- 
graded engine. 

Compact Propulsion System Model Estimates 

The CPSM estimates and associated measurements 
will be compared for the end-to-end estimation pro- 
cess. Flight condition and measurement bias effects on 
the CPSM estimates will also be reviewed. The pa- 
rameters to be compared are N\C2, N 2, AJ, PI A, 
FT IT, TT2.5, PT2.5, and 7T3. The measurements 
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for FT A , 7T2.5, and TT3 are based on single-point 
sensors, whereas the CPSM estimation produces av- 
eraged values for these parameters. A comparison of 
the CPSM estimate for net thrust is made with thrust 
stand data. 

The results of the KF estimation are an intermedi- 
ate step to the objective of the PSC algorithm. Suc- 
cessful PSC operation requires accurate CPSM esti- 
mates for such crucial parameters as SM F, FT IT , and 
FNP. Comparison of measured to estimated FT IT 
will yield an assessment of model accuracy. As for 
the SMF and FNP , no associated in-flight measure- 
ments were made and as such, these parameters must 
be estimated; this was done with established thermo- 
dynamic relationships. Inputs to the SMF and FNP 
calculations are based on measurements rather than 
estimates, when the measurements arc available. By 
comparing these measurements with the associated es- 
timates, overall model accuracy may be further sur- 
mised. However, absolute model accuracy of the SM F 
and FNP estimates cannot be realized without ad- 
ditional instrumentation and knowledge about the fi- 
delity of the thermodynamic relations. 

Baseline Results 

Data gathered on the degraded engine at a MIL- 
power setting for an acceleration at the design altitude 
of 30,000 ft serve as baseline for comparison. The de- 
graded engine was chosen for this analysis because it 
had more instrumentation than the refurbished engine. 
Although the results are from testing of a degraded en- 
gine, the CPSM should produce results independently 
of engine deterioration. The CPSM-estimated values 
for the five measured inputs N1C2, A r 2, AJ, FT A , 
FTH\ and estimated values for PT2.5, TT2.5, and 
TT3 are shown in Fig. 10 along with the correspond- 
ing measured values. Compared also are the CPSM- 
estimated fan airflow W ACC and the independent 
DEEC calculation. The following paragraphs discuss 
the CPSM estimates of Fig. 10: 

Fan Rotor Speed, Corrected to Station 2. The 

decreasing trend of N\C2 compares favorably with the 
measured value over the acceleration. Near 0.88 Mach, 
the slope of measured N\C2 decreases and the estimate 
does not change slope until approximately Mach 0.95, 
indicating some apparent lag. The estimated absolute 
value of N\C2 is as much as 1.5-percent greater than 
the measurement. 

Compressor Rotor Speed. The trend of the esti- 
mated N 2 does not agree as well as expected with its 
measurement, although neither vary significantly with 
Mach number. The difference in absolute value for N 2 
between the estimated and measured value is less than 
1 percent. 


Nozzle Throat Area. The estimated AJ is up to 
2.5-pcrcent less than the measured value. The esti- 
mated AJ is actually an effective nozzle area, related 
to the physical AJ by the discharge coefficient which 
can be quite large. As such, it can be expected that 
the effective nozzle area may be up to 10-percent less 
than the physical area. Since the estimate is less than 
the measurement, the trend is in the right direction. 

Total Pressure at the High-Pressure Turbine 
Inlet. Trends of the estimated PTA agree well with 
the measurement. The estimated value for PTA is 
always low and as much as 4-percent less than the 
measured value. 

Fan Turbine Inlet Temperature. The CPSM 
estimate of the FT IT agrees to within 20 °R of the 
measurement. Overall the trends agree well, but the 
FTIT estimate appears to lag the measurement. 

Total Pressure at the Compressor Inlet. 

Trends of the estimated PT2.5 agree very well with 
the measurement, although the value is low over the 
Mach range. The estimated PT2.5 is never more than 
2-percent less than the measured value. 

Total Temperature at the Compressor Inlet. 

The estimated 7T2.5 does a fair job in tracking the 
dynamics seen in the measurement. At the lower Mach 
numbers, the estimated 7T2.5 is 2.5-percent less than 
the measured value. 

Compressor Discharge Total Temperature. 

The estimated 77 '3 agrees well with the trends of the 
measurement. A relatively constant offset error of no 
more than 40° or 3 percent is observed between the 
estimated and measured values of TT 3. 

Digital Electronic Engine Control Calculated 
Airflow. The estimated W ACC is up to 2- percent 
greater than the DEEC calculated airflow 7 . The DEEC 
airflow is limited to a maximum of 246 Ib/sec for most 
of the acceleration, but the CPSM estimates airflow up 
to 4.5 lb/sec higher midway through the acceleration. 
Near the end of the acceleration, both airflow estimates 
are reduced and in better agreement. 

The majority of the CPSM estimates lie within a 
moderately accurate 3-percent error band over the high 
subsonic Mach region at 30,000 ft. Discrepancies in the 
estimates might be the result of either a lack of engine 
modeling fidelity, including installation and instrumen- 
tation errors, or unmodeled measurement biases. 

Measurement Bias Results 

Data were gathered on the refurbished engine at a 
MIL-power setting for an acceleration at 30,000 ft for 
two different bias configurations. The same bias sets 
used for analyzing the KF estimation are used here 
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(the final bias set 1: +70° FT IT, +2.0 in 2 AJ , and 
+ 180 lb/hr WF measurement biases; and the prelim- 
inary bias set 2: +40° FT IT and -2.9 in 2 AJ mea- 
surement biases). The C PS M- estimated values for the 
five measured inputs N\C2, N2, AJ, PTA , FTIT are 
shown in Fig. 11 along with the corresponding mea- 
sured values for both bias configurations. The CPSM- 
estimated airflow is also overplotted with the produc- 
tion DEEC-calculated airflow. 

The magnitude of estimation error for the N\C2 
seems to be more accentuated with the final bias set 1, 
whereas the preliminary bias set 2 displays less consis- 
tent behavior in the higher Mach region. The N\C2 
estimate for bias set 1 is as much as 185 rpm less than 
the measurement. The magnitude of estimation error 
for the AJ with the preliminary bias set 2 is almost 
twice as great as that with the final bias set 1. Bias 
set 2 results in an estimate as much as 10 in 2 less than 
the measured value. The magnitude of estimation er- 
ror for the FTIT is slightly more with bias set 2 than 
bias set 1. The difference of the CPSM-estimated to- 
tal engine airflow, W ACC, and the DEEC-calculated 
W ACC is greatest with bias set 2. 

Overall, the final bias set 1 produces more accurate 
CPSM estimates. In particular the estimates for AJ, 
FTIT , and W ACC appear to benefit from bias set 1. 
The magnitude of estimation error for N 2 and the PTA 
is about the same for either bias set. 

Altitude Results 

The effect of altitude on the CPSM estimation is an- 
alyzed with data gathered on the degraded engine at a 
MIL-power setting for accelerations at 15,000, 30,000, 
and 45,000 ft. The CPSM-estimated values for the five 
measured inputs N\C2, N 2, AJ, PTA, FTIT, and es- 
timated values for PT2.5, TT2.5, TT3, and W ACC 
are shown in Figs. 12 and 13 for the accelerations at 
15,000 and 45,000 ft, respectively, and from Fig. 10, 
shown previously, at 30,000 ft along with the corre- 
sponding measured values. 

The magnitude of estimation error for both the 
N \C2 and the PTA decreases with increasing altitude. 
The amount of PTA estimation error ranges from 
3 percent at 45,000 ft to 5.5 percent at 15,000 ft. The 
magnitude of estimation error for the AJ also decreases 
with increasing altitude. 

Very good estimates arc made for the N 2 and the 
FTIT. The maximum estimation error observed was 
less than 1 percent for both estimates, and is minimized 
to nearly zero at 15,000 ft. 


The errors for the estimated PT2.5, TT2.5, and TT 3 
display no clear trends with altitude. The estimation 
error for TT3 is about 3 percent throughout the three 
altitudes. 

The differences in the CPSM and DEEC total cor- 
rected engine airflow, W ACC, decrease with increasing 
altitude. At 45,000 ft, the estimated W ACC is approx- 
imately 1-percent greater than the DEEC airflow and 
about 2-percent greater at 15,000 ft and 30,000 ft. 

Good CPSM estimates are found for the MIL-power 
setting at all the flight-tested altitudes. Unexpectedly, 
the most accurate estimation occurs at altitudes other 
than the design condition, 30,000 ft. Estimation ac- 
curacy is defined as the largest observed difference be- 
tween each estimate and its corresponding sensor read- 
ing, at the same time noting the sensor accuracy. The 
table below presents, for each estimated parameter, the 
largest observed difference between measurement and 
estimate along with the sensor accuracy. The differ- 
ences represent the combined effects of sensor inaccu- 
racies and modeling uncertainty. Fortunately, the least 
accurate estimate, PTA, is also measured, and the mea- 
surement is used by the CPSM for other estimation 
outputs. The error for AJ estimate includes the dif- 
ferences between the geometric measurement and the 
model-calculated effective area. The TT3 estimation 
error did not vary with altitude; this probably repre- 
sents a discrepancy due to single- point sensor measure- 
ment not representing the average value of the engine. 
The PT2.5 and 77 2.5 estimation error probably con- 
tains effects from the assumption in the model of stan- 
dard day atmospheric conditions. 


Approximate CPSM estimation accuracy. 


Estimated 

parameter 

Sensor 

accuracy, percent 

Difference, 

percent 

N\C2 

±0.5 

1.0 

N2 

±0.5 

-1.0 

AJ 

±1.5 

-2.0 

PTA 

±3.0 

-4.0 

FTIT 

±0.5 

-0.5 

FT2.5 

±5.0 

-2.0 

TT2.5 

±0.5 

1.0 

7T3 

±0.5 

3.5 


Thrust Stand Results 


Figure 14 shows the results from a thrust stand run 
at a MIL PLA setting for the degraded engine. The 
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estimated thrust is approximately 4-percent higher 
than measured. The difference may be attributable to 
error in other PSC estimates that are used as input to 
the thrust calculation, or error in the modeled thrust 
equation. However, it is more critical for the PSC al- 
gorithm to closely match trends in thrust rather than 
the absolute values. As seen in the figure, the CPSM 
estimation performs the tracking task very well, but 
with an apparent lag of about 3 sec. Additional thrust 
stand results are presented in Ref. 7. 

Sensitivity of Model to Kalman-Filter Inputs 

Results were presented in Fig. 6 for the KF esti- 
mates produced in flight for the refurbished engine for 
an OPR of 22 (60° PLA setting) in a cruise at Mach 0.9 
and 30,000 ft. It was noted previously that the engine 
component deviation parameters were expected to be 
zero, but in real time they were not. At the 60° PLA 
setting, the KF estimates were relatively small com- 
pared with the other power settings, but the effect of 
these small numbers on the follow-on CPSM estima- 
tion was unclear. To obtain a feeling for the impor- 
tance of the engine component deviation parameters, 
the onboard CPSM model was executed postflight with 
zeroes used as the KF input to the CPSM. The result- 
ing postflight CPSM outputs are compared with the 
real-time CPSM outputs and associated measurements 
in Fig. 15. 

The estimates produced with and without the real- 
time KF values for the measured inputs N\C2, N 2, 
AJ, PTA , FT IT are overplotted with the correspond- 
ing measurements. The estimated outputs SMF and 
FNP arc also overplotted for the two cases. In 
addition, the CPSM airflow W ACC is shown for both 
cases plotted with the D EEC- calculated W ACC. 

All estimates except FT IT increase when zeroes 
are used as input from the KF to the CPSM. The 
estimated nozzle area, zl./, increases by as much as 
5 percent and the estimated FT IT decreases by about 
3 percent without the real-time KF input. The real- 
time estimates, in general, are in good agreement with 
the measurements. The estimated airflow is about 
7-percent greater with zeroes as KF input. 

The effect of the KF input on the resulting SMF 
calculation is sizable and less conservative; without the 
real-time inputs, the SMF is reduced by about 50 per- 
cent. At least some of this decrease is caused by the 
increase in estimated airflow. With zeroes substituted 
for the real-time engine component deviation parame- 
ters, the resulting thrust calculation is about 3-percent 
greater. 


For the flight test point considered, all real-time KF 
deviation parameters were negative suggesting the re- 
furbished engine is less efficient than the engine model. 
In fact, the effect on the estimated engine parameters 
agrees with this hypothesis. Estimated FNP for the 
engine model is 3 percent greater than estimated for 
the refurbished engine. In addition, the engine model, 
with zero component deviations, operates at a lower es- 
timated FTIT than for the refurbished engine. Clearly, 
even seemingly small KF estimates play an important 
role in the PSC estimation process. 

Concluding Remarks 

Test results show the propulsion system parameter 
estimation process to be successfully operating and 
producing reasonable engine estimates over the entire 
steady-state subsonic flight envelope. Variations with 
flight condition, engine degradation, and measurement 
bias indicate the sensitivity of the performance seeking 
control algorithm and estimation process. 

Comparisons with flight measurements indicate that 
the estimated inputs to the Kalman filter and compact 
propulsion system model are reasonable. Well-behaved 
estimates are produced that accurately reflect the state 
of the engine. The airflow estimation produces higher 
maximum values than those predicted by the digital 
electronic engine control calculation. Thrust stand re- 
sults show good correlation between the performance 
seeking control estimated and measured thrust. 

Reynolds effects, hardware discrepancies, and engine 
dynamics all contribute to performance seeking control 
modeling error. The sensitivity of the Kalman-filter 
engine component deviation parameters to unmodeled 
altitude and measurement bias effects indicate the need 
for improved modeling techniques. The Kalman-filter 
component deviation parameters do not accurately re- 
flect known engine degradation. The propulsion system 
model estimation is quite sensitive to small Kalman- 
filter estimates. 
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Fig. 10. Compact propulsion system model estimates compared with degraded engine measurements for a military- 
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Fig. 11. Compact propulsion system model estimates compared with measurements for 
with different biases during a military-power acceleration at 30,000 ft. 
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